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Abstract.  The  IAU  Working  Group  on  Precession  and  the  Equinox  looked  at  sev¬ 
eral  solutions  for  replacing  the  precession  part  of  the  IAU  2000A  precession-nutation 
model,  which  is  not  consistent  with  dynamical  theory.  These  comparisons  show  that  the 
(Capitaine  et  al.,  Astron.  Astrophys.,  412,  2003a)  precession  theory,  P03,  is  both  consis¬ 
tent  with  dynamical  theory  and  the  solution  most  compatible  with  the  IAU  2000A  nuta¬ 
tion  model.  Thus,  the  working  group  recommends  the  adoption  of  the  P03  precession 
theory  for  use  with  the  IAU  2000A  nutation.  The  two  greatest  sources  of  uncertainty 
in  the  precession  theory  are  the  rate  of  change  of  the  Earth's  dynamical  flattening,  A  J2, 
and  the  precession  rates  (i.e.  the  constants  of  integration  used  in  deriving  the  precession). 
The  combined  uncertainties  limit  the  accuracy  in  the  precession  theory  to  approximately 
2  mas  cent-2. 

Given  that  there  are  difficulties  with  the  traditional  angles  used  to  parameterize  the 
precession,  za,  Ca,  and  9a,  the  working  group  has  decided  that  the  choice  of  parame¬ 
ters  should  be  left  to  the  user.  We  provide  a  consistent  set  of  parameters  that  may  be 
used  with  either  the  traditional  rotation  matrix,  or  those  rotation  matrices  described  in 
(Capitaine  et  al.,  Astron.  Astrophys.,  412,  2003a)  and  (Fukushima  Astron.  J.,  126,  2003). 

We  recommend  that  the  ecliptic  pole  be  explicitly  defined  by  the  mean  orbital  angular 
momentum  vector  of  the  Earth-Moon  barycenter  in  the  Barycentric  Celestial  Reference 
System  (BCRS),  and  explicitly  state  that  this  definition  is  being  used  to  avoid  confusion 
with  previous  definitions  of  the  ecliptic. 
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Finally,  we  recommend  that  the  terms  precession  of  the  equator  and  precession  of  the 
ecliptic  replace  the  terms  lunisolar  precession  and  planetary  precession,  respectively. 

Key  words:  precession  and  the  ecliptic,  reference  systems 


1.  Introduction 

Precession  or,  more  precisely,  precession  of  the  equinox  is  the  result  of 
the  motions  of  two  planes.  The  first  plane  is  that  of  the  Earth’s  equa¬ 
tor.  The  second  is  the  ecliptic,  the  mean  plane  of  the  orbit  of  the  Earth- 
Moon  barycenter1  about  the  Sun.  These  two  planes  have  been  chosen 
because  the  equinox  has  historically  provided  a  convenient  fiducial  point 
in  the  observation  of  the  heavens  and  the  passage  of  time.  For  exam¬ 
ple,  the  calendar  year  is  tuned  to  follow  the  tropical  year  from  equinox 
to  equinox  rather  than  another  definition  of  the  year  such  as  perihelion 
passage  or  the  revolution  of  the  Earth  about  the  Sun  in  inertial  space. 
These  planes  are  also  both  dynamically  involved  in  the  motion  of  the 
Earth’s  pole. 

For  the  purposes  of  this  report,  an  inertial  space  or  an  inertial  coordinate 
system  is  one  in  which  the  space  coordinate  grids  show  no  global  rotation 
with  respect  to  a  set  of  distant  extragalactic  objects.  This  report  is  specifi¬ 
cally  concerned  with  two  coordinate  systems.  The  first  coordinate  system 
is  the  Barycentric  Celestial  Reference  System  (BCRS)  with  its  origin  at  the 
solar  system  barycenter  and  its  axes  aligned  oriented  to  match  the  Inter¬ 
national  Celestial  Reference  System  (ICRS).  The  second  coordinate  system 
is  the  Geocentric  Celestial  Reference  System  (GCRS)  with  its  origin  at  the 
center  of  mass  of  the  Earth. 

In  the  past,  the  motion  of  the  Earth’s  equator  has  been  called  lunisolar 
precession,  while  the  motion  of  the  ecliptic  has  been  called  planetary  preces¬ 
sion.  The  names  of  these  components  are  based  on  the  dominant  forces  for 
each  of  these  motions.  However,  the  accuracy  with  which  the  precession  can 
be  measured  has  reached  the  point  where  the  contribution  of  the  planets  to 
the  motion  of  the  Earth’s  equator  is  significant.  Thus,  the  name  lunisolar 
precession  has  become  misleading.  Fukushima  (2003)  proposed  renaming 
lunisolar  precession  equator  precession  and  planetary  precession  ecliptic  pre¬ 
cession  to  describe  more  accurately  these  two  components  of  precession. 


'The  ecliptic  is  defined  as  the  mean  plane  of  the  Earth’s  orbit  around  the  Sun  (The 
Astronomical  Almanac  2006).  However,  for  reasons  described  in  Hilton  (2006)  it  is  pref¬ 
erable  to  use  the  mean  orbital  angular  momentum  vector  of  the  Earth-Moon  barycenter. 
This  change  of  definition  is  included  in  recommendation  5  of  this  report. 
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Similarly,  Capitaine  et  al.  (2003a)  proposed  the  terms  precession  of  the  equa¬ 
tor  and  precession  of  the  ecliptic.  Considering  that  the  term  equator  preces¬ 
sion  does  not  differentiate  between  motion  of  the  equator  and  motion  along 
the  equator,  this  report  will  adopt  the  terms  precession  of  the  equator  and  pre¬ 
cession  of  the  ecliptic.  Further,  we  recommend  that  these  terms  be  adopted 
for  general  use. 

Since  its  adoption,  it  has  become  apparent  that  the  IAU  1976  theory 
of  general  precession  (Lieske  et  al.,  1977,  henceforth  Lieske)  is  in  error 
by  approximately  300  mas  cent-1,  where  lmas  =  O/'OO  1  and  the  century 
(cent2)  consists  of  36,525  Julian  days  Terrestrial  Time  (TT).  Williams  (1994) 
showed  that  in  addition  to  the  precession  in  longitude  there  should  also  be 
a  secular  motion  in  the  obliquity  of  the  Earth  which  he  estimated  to  be 
about  —24  mas  cent-1.  This  motion  in  latitude  is  caused  by  the  slight  incli¬ 
nation  of  the  lunar  orbit  to  the  ecliptic  when  averaged  over  the  period  of 
its  node.  More  recent  estimates  of  the  motion  in  obliquity  are:  Bretagnon 
et  al.  (2003),  —25  mas  cent-1;  Capitaine  et  al.  (2003a),  —26  mas  cent-1;  and 
Fukushima  (2003),  —23  mas  cent-1. 

The  precession  and  nutation  of  the  Earth  are  most  accurately  observed 
using  Very  Long  Baseline  Interferometer  (VLBI)  observations.  These  obser¬ 
vations  are  only  sensitive  to  the  linear  portion  of  the  precession  and  insen¬ 
sitive  to  the  ecliptic.  Thus,  the  higher  order  coefficients  of  the  precession 
theory  along  with  the  orientation  of  the  ecliptic  must  be  derived  from 
dynamical  theory.  When  the  IAU  2000  precession-nutation  theory  (IERS, 
2004)  was  adopted  (IAU,  2001)  the  emphasis  of  the  analysis  was  on  the 
periodic  nutations  and  correcting  the  linear  portion  of  the  precession  VLBI 
observations.  The  effect  of  these  changes  on  the  higher-order  terms  in 
the  precession  theory  was  ignored.  Ignoring  the  higher-order  terms  results 
in  an  error  in  the  precession  of  about  6.4  mas  cent-2  in  longitude  and 
0.01  mas  cent-2  in  obliquity.  Thus,  the  precession  theory  was  not  consistent 
with  dynamical  theory. 

Fukushima  (2003)  showed  that  the  values  of  £a  and  za,  two  of  the  tra¬ 
ditional  angles  for  parameterizing  the  precession,  are  complementary  and 
highly  dependent  on  the  precise  values  that  are  adopted  for  the  offset 
between  the  dynamical  frame  and  the  GCRS  at  J2000.0.  Thus,  they  are 
unsuited  to  polynomial  representation. 

The  ecliptic  in  use  was  defined  by  Lieske  using  a  method  similar  to  that 
of  Newcomb  (1894).  These  yield  an  ecliptic  which  follows  the  geometrical 

2The  IAU  style  book  (Wilkins,  1987)  discourages  the  use  of  century  as  a  unit;  how¬ 
ever,  it  is  frequently  used  for  slowly  varying  angles  in  celestial  mechanics  applications.  Nor 
is  there  a  consensus  for  its  abbreviation.  The  abbreviation  used  here,  cent,  was  chosen  as 
the  one  found  most  often  in  a  sample  of  both  British  and  American  English  dictionaries. 
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path  of  the  Earth-Moon  barycenter.  However,  more  recent  work  has 
used  instead  an  ecliptic  that  follows  the  Earth-Moon  barycenter’s  orbital 
angular-momentum  vector  in  the  BCRS,  yielding  a  slightly  different  result 
due  to  the  small  contribution  from  the  rotation  of  the  ecliptic  itself  that 
is  included  in  the  former  method  (see  Standish,  1981).  Thus  both  the 
geometrical,  often  referred  to  as  “rotating”,  and  inertial  definitions  of  the 
ecliptic  have  been  used,  creating  confusion. 

The  IAU  Working  Group  on  Precession  and  the  Ecliptic  was  formed  at 
the  XXYth  General  Assembly  of  the  IAU  in  Sydney,  Australia  to  address 
these  topics  and  make  recommendations  regarding  them  to  the  IAU.  This 
report  constitutes  the  findings  of  that  working  group. 


2.  Choice  of  Precession  Theory 

Four  recent,  high  precision  theories  for  the  precession  and  the  definition 
of  the  ecliptic  (Bretagnon  et  al.,  2003;  Capitaine  et  al.,  2003a;  Fukushima, 
2003;  Harada  and  Fukushima,  2004)  were  considered.  After  examina¬ 
tion  of  the  four  theories,  the  working  group  agreed  that  precession  of 
the  equator  in  the  Capitaine  et  al.  (2003a)  precession  theory,  P03,  is 
consistent  with  dynamical  theory  and  recommends  it  as  the  high  preci¬ 
sion  precession  theory  most  suited  for  use  with  the  IAU  2000A  nuta¬ 
tion.  The  details  of  the  comparison  of  the  four  theories  can  be  found 
in  Capitaine  et  al.  (2004).  The  two  greatest  sources  of  uncertainty  in 
the  precession  theory  are  the  rate  of  change  of  the  Earth’s  dynamical 
flattening,  A  J2,  and  the  precession  rates.  The  portion  of  the  uncertainty 
from  the  uncertainty  in  A/2  limits  the  accuracy  in  the  precession  theory  to 
approximately  1.5  mas  cent-2  (Bourda  and  Capitaine  2004).  And  the  uncer¬ 
tainty  in  the  precession  rates  is  approximately  the  same.  Thus,  the  com¬ 
bined  uncertainties  limit  the  accuracy  in  the  precession  theory  to  about 
2  mas  cent-2. 

While  it  serves  a  useful  fiducial  purpose,  dynamical  uses  for  the  mod¬ 
ern  ecliptic  are  limited.  Thus,  certain  arbitrary  decisions  required  to  realize 
the  ecliptic  do  not  impair  it.  The  details  of  these  considerations  are  given 
in  Hilton  (2006).  Thus,  the  working  group  recommends  the  precession 
of  the  ecliptic  included  in  the  P03  precession  theory.  More  specifically, 
we  recommend  that  the  ecliptic  pole  should  be  explicitly  defined  by  the 
mean  orbital  angular  momentum  vector  of  the  Earth-Moon  barycenter 
in  the  BCRS  to  simplify  the  dynamics.  We  also  recommend  that  both 
the  definition  used  and  the  process  by  which  the  ecliptic  has  been  deter¬ 
mined  be  made  explicit  when  any  future  definition  is  adopted,  to  avoid 
confusion. 
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3.  Parameterization 

There  is  no  unique  method  of  parameterizing  the  precession.  Williams 
(1994),  for  example,  gives  three  methods  to  rotate  from  the  equator  and 
equinox  of  J2000.0  to  the  equator  and  equinox  of  date  in  addition  to  the 
method  used  by  Lieske  (1979).  The  particular  parameterization  desired  may 
depend  on  extraneous  details  such  as  whether  old  computer  code  needs  to 
be  maintained  or  replaced.  Thus,  we  have  chosen  to  provide  the  data  nec¬ 
essary  for  four  different  parameterizations: 

•  The  traditional  parameterization  used  by  Newcomb  and  Lieske. 

•  The  parameterization  recommended  by  Capitaine  et  al.  (2003a)  which 
cleanly  separates  precession  of  the  equator  from  precession  of  the  eclip¬ 
tic. 

•  The  parameterization  developed  by  Fukushima  (2003)  which  allows  for 
flexible  switching  between  the  mean  equator  and  equinox  of  J2000.0  and 
the  GCRS  as  its  initial  reference  system. 

•  The  parameterization  to  compute  the  precession-nutation  in  the  new 
paradigm  adopted  by  the  IAU  with  respect  to  the  Celestial  Intermediate 
Origin  (CIO)3. 

The  first  two  parameterizations  use  the  position  of  the  mean  equator  and 
equinox  of  J2000.04  as  their  initial  reference  frame.  However,  many  applica¬ 
tions  nowadays  use  the  GCRS  as  the  initial  reference  frame.  To  start  from  the 
GCRS,  the  coordinates  used  in  the  first  two  parameterizations  must  first  be 
rotated  using  a  constant  bias  matrix  (top  line  of  Expression  4  in  Capitaine 
et  al.,  2003a).  The  third  parameterization  may  start  from  either  the  mean 
equator  and  equinox  of  J2000.0  or  the  GCRS.  The  description  provides  the 
details  on  how  the  parameters  change  between  the  mean  equator  and  equi¬ 
nox  of  J2000.0  and  the  GCRS.  The  fourth  parameterization  is  meant  only 
to  be  used  with  the  GCRS  as  the  initial  reference  frame. 

The  polynomial  coefficients  for  all  the  angles  for  these  four  parameter¬ 
izations  derived  from  the  P03  precession  are  given  in  Table  I.  The  first  15 
expressions  are  from  Capitaine  et  al.  (2003a).  The  remaining  six  expres¬ 
sions  are  the  result  of  applying  the  P03  precession  to  the  angles  defined  by 
Fukushima  (2003). 


3The  IAU  Division  I  Working  Group  on  Nomenclature  for  Fundamental  Astronomy 
is  recommending  the  Celestial  Ephemeris  Origin  (CEO)  be  renamed  the  Celestial  Inter¬ 
mediate  Origin  (CIO). 

4The  position  of  the  CIP,  not  including  nutation,  at  J2000.0. 


TABLE  I 

The  polynomial  coefficients  for  the  precession  angles. 
Coefficients* 
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*Centuries  (cent.)  are  Julian  centuries  of  36,525  days  TT. 
t  The  angle  e0  =  £a  (t  =  0) . 
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3.1.  LIESKE  PARAMETERIZATION 

The  Lieske  (1979)  precession  matrix  is: 

R(za,#a,?a) 

=  ^3  (—90°  —  za)Ri  (9a)  R3  (90°  —  fA) 

=  R3{-za)R2(9a)R}(-U) 

cosza  cos$a  cos£a — sinzA  sin  £A  —  coszACOS$AsmfA— sinzAcosfA  — coszAsin(?A 
sinzACOS@ACOSfA+cosjAsmfA  —  sinzAcos0AsinfA+coszACOSfA  — sinzAsin6>A 
sin^ACOSfA  — sin^AsinfA  cos0A 

(1) 

where  R\,  Ri,  and  /p  are  the  basic  rotation  matrices  around  the  axes  of  a 
right  handed  coordinate  system, 

(10  0  \  / cost  0  — sinr\  /  cost  sinT  0\ 

0  cost  sinT  ,  R2(r)=  0  1  0  ,  and  R3(t)=  — sinT  cost  0 

0  —sinT  cost  /  \  sinT  0  cost  )  \  0  0  1/ 

(2) 

The  angles,  shown  in  Figure  1,  parameterizing  the  rotation  are:  6\,  the  arc 
connecting  the  mean  pole  of  J2000.0,  Po,  with  the  mean  pole  of  date,  P,: 
£a,  the  angle  with  its  vertex  at  Po  from  the  equinoctial  colure  of  J2000.0  to 
6»a;  and  za,  the  angle  with  its  vertex  at  P,  from  dA  to  the  equinoctial  colure 
of  date.  This  might  be  considered  the  traditional  way  of  parameterizing  the 
precession  and  is  imbedded  in  many  software  applications  that  implement 
the  precession.  The  angles  £a,  Pa-  and  za  can  be  determined  from  coa,  the 
arc  from  the  ecliptic  pole  of  J2000.0  to  the  mean  pole  of  date;  ipA,  the 
angle  between  eo,  the  arc  from  the  ecliptic  pole  of  J2000.0  to  the  mean  pole 
of  J2000.0,5  and  coa,  tza,  the  obliquity  of  the  ecliptic  of  date  on  the  ecliptic 
of  J2000.0;  and  nA,  the  longitude  of  the  ascending  node  of  the  ecliptic  of 
date  on  the  ecliptic  of  J2000.0  (see  Figure  1).  The  last  two  quantities,  tta 
and  nA,  are  given  in  terms  of  the  canonical  variables  PA  =  sin  tta  sin  11 A 
and  <2a  —  skuTACOS  FIa-  The  relations  between  the  known  quantities  and 
?a,  0A,  and  zA  are 

cos  9a  =  cos  eo  cos  co  a  +  sin  eo  sin  ooa  cos  x/za  , 
cosCA  =  sin&)Asini/'A/sin0A,  and  (3) 

cos(zA  -  Xa)  =  sin  €0  sin  ^A/  sin  9 a 

where  the  length  of  eo  =  coa(z  =  0),  and  xa,  the  angle  between  coA  and  6a,6  is 
sin  XA  =  sin(180°  —  FIa  —  i//A )  sin  7r/ sin  cA  (4) 

5The  angular  length  of  the  arc  e0  is  the  obliquity  of  the  ecliptic  at  J2000.0. 

6The  arc  from  the  ecliptic  pole  of  date  to  the  mean  pole  of  date.  The  angular  length  of  this 
arc  is  the  obliquity  of  the  ecliptic  of  date. 
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Figure  1.  The  angles  associating  the  mean  pole  of  J2000.0  (JD  2451545.0),  P0;  the  mean 
pole  of  date,  Pt;  the  ecliptic  pole  of  J2000.0,  Co;  the  ecliptic  pole  of  date,  C;  the  mean 
equinox  of  J2000.0,  T0;  and  the  mean  equinox  of  date,  T,  used  in  the  Lieske  and  Cap- 
itaine  et  al.  (2003)  parameterizations  for  the  precession. 

where 

coseA  =  cos&)ACOS^A  +  sin&>A  sin  7Ta  cos  (180°  -  nA  —  iAa)  (5) 


A  difficulty  arises  because  the  starting  point  for  many  applications 
is  now  the  GCRS  rather  than  mean  equator  and  equinox  of  J2000.0. 
Although  this  can  be  accommodated  in  the  traditional  angles,  Fukushima 
(2003)  and  Capitaine  et  al.,  (2003a)  have  shown  that  values  of  the  angles 
za  and  Ca  near  J2000  are  highly  sensitive,  in  a  complementary  way,  to  the 
values  adopted  for  the  offset  between  the  mean  equator  and  equinox  of 
J2000.0  and  the  GCRS  and  hence  cannot  satisfactorily  be  represented  by 
polynomials  in  t.  Thus,  unless  the  accuracy  goals  are  sufficient  to  neglect 
the  frame  bias  (about  23  mas  overall),  either  the  frame  bias  has  to  be  intro¬ 
duced  as  an  additional  rotation,  or  a  different  representation  that  does  not 
suffer  from  these  difficulties  should  be  adopted. 


3.2.  CAPITAINE  ET  AL.  PARAMETERIZATION 

Capitaine  et  al.  (2003a)  provide  a  clean  separation  between  the  precession 
of  the  equator  and  precession  of  the  ecliptic  providing  the  precession  of  the 
equator  with  respect  to  the  equator  and  equinox  of  J2000.0,  V'A  and  coa, 
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and  the  precession  of  the  ecliptic,  Pa  and  Q\,  with  respect  to  the  ecliptic 
and  equinox  of  J2000.0.  The  most  efficient  method  of  computing  the  pre¬ 
cession  matrix  using  these  parameters  is  first  to  determine  the  mean  obliq¬ 
uity  of  date,  e\,  and  xa,  the  angle  between  the  ecliptic  pole  of  J2000.0 
(epoch),  Co,  and  the  ecliptic  pole  of  date,  C,  with  its  vertex  at  P,  (see  Fig¬ 
ure  1).  The  angle  xa  is  determined  using  Equation  4  and  e\  is  determined 
using  Equation  5.  The  precession  matrix  is  then: 

/  Pn  Pi  2  Pn\ 

P(/A,  <WA,  VfA)  =  P3(XA)Pl(-«A)P3(-^A)Pl(eo)=  T21  P22  P23  (6) 

\  P31  P32  P33  / 

where  the  components  of  the  rotation  matrix  are: 

Pi  1  =  cos  xa  cos  V^a  +  sin  xa  cos  <x>a  sin  xpA , 

P12  =  (—  cos  xa  sin  xpA  +  sin  xa  cos  coa  cos  xpA)  cos  6o  +  sin  xa  sin  ooA  sin  60, 

Pi 3  =  (—  cos  xa  sin  xpA  +  sin  xa  cos  coa  cos  xpA)  sin  60  —  sin  xa  sin  a>A  cos  eo, 

P21  =  —  sin  xa  cos  xpA  +  cos  xa  cos  coa  sin  x[r A, 

P? 2  =  (sin  xa  sin  xpA  +  cos  xa  cos  coa  cos  xpA)  cos  eo  +  cos  xa  sin  coA  sin  eo,  (7) 
P2 3  =  (sin  xa  sin  x//A  +  cos  xa  cos  coa  cos  xpA)  sin  eo  —  cos  xa  sin  a>A  cos  eo, 

P31  =  sin(UAsini/fA, 

P3?  =  sin  coA  cos  x/sA  cos  €0  —  cos  coA  sin  c0 ,  and 
P33  =  sin  coA  cos  x//A  sin  60  +  cos  coA  cos  eo  ■ 


3.3.  FUKUSHIMA  PARAMETERIZATION 

Fukushima  (2003)  provides  an  alternative  parameterization  that  has  certain 
advantages  over  the  parameterization  used  by  Lieske  (1979)  or  the  separa¬ 
tion  into  precession  of  the  equator  and  precession  of  the  ecliptic  provided 
by  Capitaine  et  al.  (2003a).  This  parameterization  starts  with  P0,  and  the 
arc  to  Co  on  the  celestial  sphere  (see  Figure  2).  The  position  of  the  C  is 
then  determined  by  the  arc,  </>,  from  P0  to  C,  and  y,  the  angle  between  (j> 
and  the  arc  from  P0  to  Co  with  its  apex  at  Po.  The  position  of  the  mean 
pole  of  date  is  then  determined  from  <?a  and  ij/1 ,  the  angle  between  </>  and 
6a  with  its  apex  at  C.  The  angles  (p,  y,  and  xp  are  related  to  already  known 
angles  by: 


7The  angle  1 {r  is  not  the  same  as  the  angle  ir A.  The  former  is  measured  from  (p  to 
6a  with  the  apex  at  the  ecliptic  pole  of  date,  while  the  latter  is  measured  from  60  to  a> a 
with  the  apex  at  the  ecliptic  pole  of  J2000.0. 
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Figure  2.  The  angles  associating  the  mean  pole  of  date,  P,\  and  the  ecliptic  pole  of  date, 
C;  to  the  mean  pole  of  J2000.0,  P0;  and  ecliptic  pole  of  J2000.0,  C0  in  the  Fukushima 
(2003)  parameterization  for  the  precession. 

cos  4>  —  cos  7Ta  cos  —  sin  jt\  sin  <?0  cos  nA  (8) 


sin  7ta  sin  11a  Pa 

sin  y  = - ; - =  — — 

sm  cp  sin  cp 


(9) 


sin  6 a  cos  za 

sin  t fr  = - ; - 

sm  </> 

The  precession  matrix  is  then: 


(10) 


P(e  a,  iM,  y)  =  Ri(-eA)R3(-^)Ri((p)R3(y)  = 


P\\ 

Pl2 

P 1 3 

P21 

P22 

^23 

P31 

P32 

T33 

(11) 


where  the  components  of  the  rotation  matrix  are: 


Pu  =  cos  t/r  cosy  +  sint/rcos(/>siny, 

P\2  =  cos  t/r  sin  y  —  sin  x/r  cos  4>  cos  y , 

Pi 3  =  —  sin  sirup, 

P21  =  cos  e\  sin  1 (r  cos  y  —  (cos  cos  1 jr  cos  4>  +  sin  sin  4>)  sin  y, 

P22  =  cosca  sint/r  sin  y  +  (cos cos  xp  cos 4>  +  sin sin 4>)  cosy,  (12) 
P23  =  cos  eA  cos  t/r  sin  cp  —  sin  cos  <p, 

P31  =  sin  6a  sin  t/r  cos  y  —  (sin  6A  cos  x[r  cos  cp  —  cos  6A  sin  cp)  sin  y, 

P32  =  sin  6a  sin  i/r  sin  y  +  (sin  eA  cos  i/r  cos  </>  —  cos  e\  sin  cp)  cos  y ,  and 
P33  =  sin  6a  cos  t/r  sin  <p  +  cos  6A  cos  cp . 
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For  some  applications,  the  initial  reference  frame  is  the  GCRS  rather 
than  the  mean  equator  and  equinox  of  J2000.0.  Using  the  GCRS  as 
the  initial  reference  frame  requires  three  changes.  The  relevant  angles  are 
shown  in  Figure  3. 

First,  in  Equation  8,  the  arc  eo  is  replaced  by  eo  +  <5c0  where  <5eo  is  the 
difference  between  the  mean  obliquity  of  the  ecliptic  at  J2000.0  and  the  arc 
length  from  Co  to  the  pole  of  the  GCRS,  Rgcrs-  Thus, 

cos0gcrs  =  cos7TaCOs(6o  +  <$eo)  —  shuta  sin(eo  +  <5eo)  cos  nA.  (13) 

Second,  unlike  the  mean  equator  and  equinox  of  J2000.0,  the  node  of 
the  ecliptic  of  J2000.0  on  the  GCRS  equator  is  offset  0."052928  in  right 
ascension.  Thus,  an  additional  initial  rotation,  = —0. "052928,  about  the 
z-axis  is  required  to  first  align  the  origin  of  the  GCRS  with  the  node  of 
the  ecliptic  on  the  equator  of  the  GCRS,  and  the  equation  for  y  becomes 


(14) 


Third,  let  8%[/  be  the  angle  between  the  arc  to  P,  and  the  arc  to  Rgcrs 
with  its  vertex  at  C,  and  8x//0  be  the  angle  between  the  arc  to  Po  and  the 
arc  to  P gcrs  with  its  vertex  at  Co.  The  angles  <5  and  Sxjr  are  very  small. 
Thus,  the  cosines  of  these  angles  are  ~(l-5x  10-17).  Correct  calculation 


Co 


Figure  3.  The  angles  associating  the  mean  pole  of  date,  Pt;  and  the  ecliptic  pole  of  date, 
C;  to  the  mean  pole  of  J2000.0,  P0;  and  ecliptic  pole  of  J2000.0,  C0  in  the  Fukushima 
(2003)  parameterization  for  the  precession. 
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of  these  angles  using  a  cosine  formulation  requires  extended  precision  to 
retain  accuracy.  Instead,  it  is  preferable  to  use  the  half  angle  sine  formula, 
Equation  C3,  of  Fukushima  (2003).  The  length  of  the  arc,  S,  from  P0  to 
Tgcrs  is  then 


=  sin  60  sin(eo  +  8e  o)  sin2 


+  sin2 


(15) 


Since  the  angles  8e o  and  Sxpo  are  very  small,  the  sines  of  these  angles  may 
be  replaced  with  their  first  order  approximations  without  any  loss  of  accu¬ 
racy.  Similarly,  sin(eo  +  <5eo)  may  be  replaced  by  sin  e0  without  any  loss  of 
accuracy.  Thus, 


(16) 


The  same  half  angle  formula  is  used  to  find  8ip 

■  2  f  S\  ■  ,  ■  ,  ,  .  •  2  ,  ■  2(4*  0GCRS  \  nl 

sm  I  -  l=sm0sm(0GCRs)sin  1  —  1+ sin  I - ^ - )•  (I7) 


As  before,  the  angles  8xp  and  cp  —  cp gcrs  are  very  small.  Thus,  the  sines  of 
these  angles  may  be  replaced  with  their  first  order  approximations  with¬ 
out  any  loss  of  accuracy.  Also,  either  sin  (/>gcrs  may  be  replaced  by  sirup, 
or  sin  cp  may  be  replaced  by  sin</>GCRS  without  any  loss  of  accuracy.  Thus, 
making  the  appropriate  substitutions,  Equation  17  becomes 


2 

sin2  (p  + 


/  4>  ~  </>GCRS 

V  2 


2 


(18) 


The  value  of  <5^  is  determined  by  equating  the  right  hand  side  of  Equation 
16  with  the  right  hand  side  of  Equation  18  and  solving  for  Sip-.  The  result 
is 

-J  sin2  e08xp 2  +  8ej  -(cp-  cpocRs)2 

- .  (19) 

sirup 

The  negative  value  of  the  square  root  function  is  specifically  chosen  since 
the  rotation  about  the  ecliptic  pole  of  date  from  pole  of  the  GCRS  to 
the  Celestial  Intermediate  Pole  (CIP)  is  clockwise  as  shown  in  Figure  3. 
At  J2000.0,  cp  =  e 0  and  cp  —  </>gcrs  =  <$<u).  Tims*  8\p  =  8\p0  at  J2000.0,  and 
the  rotation  matrix  reduces  to  the  bias  matrix  between  the  GCRS  and  the 
mean  equator  and  equinox  of  J2000.0. 

Using  these  transformations,  the  pole  of  the  GCRS  may  be  substi¬ 
tuted  for  the  mean  pole  of  J2000.0  using  <5eo  =  AeGCRS  and  8x/r o  =  A^gcrs 
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(Mathews  et  al.,  2002).  Thus,  the  bias  matrix  to  rotate  from  the  GCRS 
to  the  mean  equator  and  equinox  of  J2000.0  is  incorporated  into  the  pre¬ 
cession.  The  angles  Se o  and  SxJ/o  are  small  enough  that  for  this  particular 
transformation  the  only  significant  differences  are  in  the  zeroth  and  first 
order  polynomial  coefficients  for  y,  4>,  and  x Jr . 

A  second  advantage  of  this  parameterization  is  that  the  intermediate 
frame  resulting  after  the  second  rotation,  Rx(4>),  is  with  respect  to  the 
ecliptic  of  date.  Thus,  nutation  may  be  added  simply  by  adding  the  nuta¬ 
tion  in  longitude,  Ax/r,  to  the  third  rotation  and  the  nutation  in  obliquity, 
Ae,  to  the  fourth  rotation,  that  is,  the  rotation  to  the  true  pole  of  date  is: 


NP(eA,  Ae,  xjr,  AxJ/,  0,  y)  =  Rx(-e A  -  A e)R3(-xJr  -  AxJ/)Rx(cJ>)R3(y). 

(20) 

While  this  parameterization  allows  the  precession  and  nutation  to  be 
applied  directly,  the  two  previous  parameterizations  require  nutation  to  be 
done  in  an  extra  step: 


NP  =  N(e A,  Ae,  A xJ/)P  =  Rx(-eA  -  Ae)R3(AxJ/)Rx (eA)P.  (21) 

Note  that  the  reason  that  the  nutation  matrix  has  to  be  separate  from 
the  precession  matrix  in  the  Lieske  and  the  Capitaine  et  al.  parameter¬ 
izations  is  that  the  precession  is  referred  to  the  ecliptic  pole  of  J2000.0 
while  the  nutation  is  referred  to  the  ecliptic  pole  of  date.  Were  the  nuta¬ 
tion  referred  to  the  ecliptic  pole  of  J2000.0  rather  than  the  ecliptic  pole 
of  date,  the  nutation  could  be  incorporated  into  the  Lieske  and  Capitaine 
et  al.  parameterizations  using  a  method  similar  to  that  of  the  Fukushima 
parameterization . 


3.4.  THE  CELESTIAL  INTERMEDIATE  ORIGIN  (CIO) 

The  CIO  is  intended  for  use  with  only  a  combined  bias-precession-nutation 
matrix  (Capitaine  et  al.,  2003b).  This  parameterization  is  in  terms  of  the 
two  angles  X,  in  the  x -direction  (toward  the  x-origin  of  the  GCRS,  and  Y , 
in  the  y -direction  (90°  to  the  east  of  x),  which  give  the  x  and  y-coordinates 
of  the  CIP  unit  vector  in  the  GCRS  (Figure  4).  The  third  parameter  is  5. 
Except  for  a  tiny  fixed  offset,  this  parameter  is  the  difference  between  the 
length  of  the  arcs  from  E0  to  N'  and  a  to  N'.  The  point  E0  is  the  x-origin 
of  the  GCRS,  a  is  the  CIO  of  date,  and  N'  is  the  node  of  the  equator 
of  the  CIP  on  the  equator  of  the  GCRS.  This  matrix  includes  precession, 
nutation,  coupling  between  precession  and  nutation,  and  frame  biases. 
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Figure  4.  The  Celestial  Intermediate  Origin  (CIO),  a ,  and  the  position  of  the  Celestial 
Intermediate  Pole  (CIP),  P,  with  respect  to  the  pole,  Pgcrs>  and  equator  of  the  Geocen¬ 
tric  Celestial  Reference  System  (GCRS)  and  its  origin,  S0.  The  parameter  X  is  the  com¬ 
ponent  of  the  arc  from  Pqcrs  to  P  in  the  x  direction  and  Y  is  the  component  in  the 
y  direction.  Thus  the  CIP  is  located  at  [X,  Y,  1  —  ( X 2  +  T2)1/2]  with  respect  to  the  GCRS 
located  at  [0, 0,  1].  The  third  parameter,  s,  is  the  difference  between  the  arcs  oN'  and 
£0AG 


The  bias-precession-nutation  matrix  is: 

1 -aX2  -aXY 
NPB  =  /?3(-s)-  I  -aXY  \-aY 2 


-X 

—  Y 


X 


Y  \-a(X2  +  Y2) 


where 


1  (X2  +  Y2) 

a&-  +  - - - 

2  8 

5  3 

X  =  ^  '  Xjt'  +  EE  tk(aSjkSinbk  +  aCjkCOsbk),  and 

i= 0  ./  k= 0 

5  3 

Y  =  ^  vT'  +  y,  y,  tk (dsjic  sin  fk  +  dcjk  cos  fk) 
i= 0  ./  k=0 


(22) 


(23) 


t  is  the  time  in  Julian  centuries  from  J2000.0  TT,  x,  and  y,  are  the  coeffi¬ 
cients  for  the  frame  bias  and  precession  in  X  and  Y,  respectively,  aSjk,  acjk, 
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dSjk,  and  dcjk  are  the  coefficients  for  the  nutation  and  coupling  between 
precession  and  nutation,  and  bk  and  fk  are  the  fundamental  angular  argu¬ 
ments  for  the  nutation  and  coupling  terms. 

The  values  for  the  coefficients  for  the  secular  part  of  the  relations  for  X 
and  Y  are  almost  entirely  due  to  precession  plus  the  bias  between  BCRS 
and  the  J2000  mean  pole.  These  coefficients  are  given  in  Table  I.  The 
parameter  s  is  also  represented  by  a  series,  but  it  is  more  efficient  to  derive 
s  from  the  series  for  s  +  XY/2,  also  in  Table  I.  The  equations  for  the  sec¬ 
ular  portions  of  X,  Y,  and  s  +  XY/2  are  also  given  in  Equations  49-51  of 
Capitaine  et  al.  (2003a). 

Determination  of  the  coefficients  for  the  periodic  parts  of  the  series  for 
X  and  Y  require  both  the  precession  and  the  nutation  portion  of  IAU 
2000A  (Mathews  et  al.  2002).  Thus,  developing  coefficients  for  the  periodic 
parts  of  the  series  was  deemed  beyond  the  scope  of  the  Working  Group. 
However,  only  very  slight  adjustments  are  required  to  the  IAU  2000  nuta¬ 
tion,  and  likewise  to  the  periodic  part  of  the  X  and  Y  series,  in  order  to 
become  consistent  with  P03,  and  the  periodic  part  of  s  +  XY/2  is  insen¬ 
sitive  to  this  change.  An  explanation  of  the  theory  and  the  sources  for 
machine  readable  tabulations  of  the  coefficients  compatible  with  the  IAU 
2000A  precession-nutation  model  are  on  pg.  44  of  IERS  (2004). 

The  elements  of  the  matrix  in  Equation  22  are  dimensionless.  Thus,  the 
values  of  A,  Y,  s,  and  the  coefficients  for  the  nutation  and  coupling  between 
precession  and  nutation  should  be  in  radians.  However,  these  values  are  nor¬ 
mally  given  in  arcseconds,  as  in  Table  I,  and  have  to  be  converted. 


4.  Recommendations 

The  Working  Group  on  Precession  and  the  Ecliptic,  recognizing: 

1.  the  need  for  a  precession  theory  consistent  with  dynamical  theory  and 
compatible  with  the  IAU  2000A  nutation  theory, 

2.  the  gravitational  attraction  of  the  planets  make  a  significant  contribu¬ 
tion  to  the  motion  of  the  Earth’s  equator,  making  the  terms  lunisolar 
precession  and  planetary  precession  are  misleading, 

3.  the  need  for  a  definition  of  the  ecliptic  for  both  astronomical  and  civil 
purposes,  and 

4.  recognizing  the  ecliptic  has  been  defined  both  with  respect  to  an 
observer  situated  in  inertial  space  (inertial  definition)  and  an  observer 
comoving  with  the  ecliptic  (rotating  definition), 
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recommends: 

1.  The  terms  lunisolar  precession  and  planetary  precession  be  replaced  by 
precession  of  the  equator  and  precession  of  the  ecliptic,  respectively. 

2.  The  IAU  adopt  the  P03  precession  theory,  of  Capitaine  et  al.,  (2003a, 
Astron.  Astrophys.,  412,  567-586)  for  the  precession  of  the  equator 
(Equations  37)  and  the  precession  of  the  ecliptic  (Equations  38);  the 
same  paper  provides  the  polynomial  developments  for  the  P03  primary 
angles  and  a  number  of  derived  quantities  for  use  in  both  the  equinox 
based  and  CIO  based  paradigms. 

3.  The  choice  of  precession  parameters  be  left  to  the  user. 

4.  The  recommended  polynomial  coefficients  for  a  number  of  precession 
angles  are  given  in  Table  I  of  this  report,  including  the  P03  expres¬ 
sions  set  out  in  Tables  3-5  of  Capitaine  et  al.  (2005,  Astron.  Astrophys., 
432,  355-367),  and  those  of  the  alternative  Fukushima  (2003,  Astron. 
J.,  126,  494)  parameterization;  the  corresponding  matrix  representations 
are  given  in  Equations  1,  6,  11,  and  22. 

5.  The  ecliptic  pole  should  be  explicitly  defined  by  the  mean  orbital  angu¬ 
lar  momentum  vector  of  the  Earth-Moon  barycenter  in  the  BCRS,  and 
this  definition  should  be  explicitly  stated  to  avoid  confusion  with  older 
definitions. 
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